Cumulant generating function formula of heat transfer in ballistic system with 

lead-lead coupling 



3 



Huanan Li0 Bijay Kumar Agarwalla, and Jian-Sheng Wang 
Department of Physics and Center for Computational Science and Engineering, 
National University of Singapore, Singapore 117542, Republic of Singapore 

(Dated: 24 Aug 2012) 

Based on two-time observation protocol, we consider heat transfer in a given time interval £m in 
lead-junction- lead system taking coupling between the leads into account. In view of the two-time 
observation, consistency conditions are carefully verified in our specific family of quantum histo- 
ries. Furthermore, its implication is briefly explored. Then using nonequilibrium Green's function 
(NEGF) method, we obtain an exact formula for the cumulant generating function for heat transfer 
between the two leads valid in both transient and steady-state regimes. Also, a compact formula for 
the cumulant generating function in the long-time limit is derived, for which the Gallavotti- Cohen 
(GC) fluctuation symmetry is explicitly verified. In addition, a brief discussion of Di Ventra's trick 
regarding whether the effect of the repartitioning procedure of the total Hamiltonian on nonequi- 
libium steady-state current fluctuation exists is given. All kinds of properties of nonequilibrium 
current fluctuations such as the fluctuation theorem in different time regimes could be readily given 
according to these exact formulas. 
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I. INTRODUCTION 

The physics of nonequilibrium many-body systems is 
one of the most rapidly expanding areas of theoretical 
physics. In the combined field of non-equilibrium state 
and statistics, the distribution of transferred charges in 
electronic case or heat in phononic case, the so-called 
full counting statistics (FCS), plays an important role, 
according to which we could understand the general fea- 
tures of currents and their fluctuations. Also, it is well 
known that the noise generated by nanodevices contains 
valuable information on microscopic transport processes 
not available from only transient or steady current. In 
FCS, the key object we need to study is the cumulant 
generating function (CGF) which presents high order cor- 
relation information of the corresponding system for the 
transferred quantity. 

The study of the FCS started from the field of elec- 
tronic transport pioneered by Levitov and Lesovik, who 
presented an analytical result for the CGF in the long- 
time limiti. After that, many works followed in electronic 
FCS^r— , while much less attention is given to phonon 
transport. Saito and Dhar are the first ones to bor- 
row this concept to thermal transport^. Later, Ren and 
his co-workers gave results for two-level systems^. And 
very recently, transient behavior and long-time limit of 
CGF have been obtained in lead-junction-lead harmonic 
networks both classically and quantum-mechanically us- 
ing Langevin equation method and NEGF method, 
respectively^— . Experimentally the FCS in electronic 
case has been carefully studied and the cumulants to very 
high order have been successfully measured in quantum- 
dot system o 11 ' 12 . In principle, similar measurements 
could be carried out for thermal transport, e.g., in a 
nanoresonator system. 

Through modern nanoscalc technology, small junction 



is easily realized such as in certain nanoscalc systems, 
for instance, a single molecule or, in general, a small 
cluster of atoms between two bulk electrodes. In that 
case, the electrode surfaces of the bulk conductors may 
be separated by just a few angstroms so that some finite 
electronic coupling between the two surfaces is present 
taking into account the long-range interaction. In or- 
der to study steady current in this situation, Di Ventra 
suggest that we can choose our "sample" region (junc- 
tion) to extend several atomic layers inside the bulk elec- 
trodes where screening is essentially complete so that the 
above coupling could be negligible^. This trick has been 
checked in some restricted case in a recent worki^. How- 
ever, whether this intuitively reasonable trick applies to 
all the higher cumulants of heat transfer in steady state 
is still a question. Therefore we will briefly discuss this 
question at the end of the article. And obviously this 
trick can not be applied to study transient behavior of 
all the cumulants of heat transfer. Thus, in this work, 
we want to construct the CGF of heat transfer in a gen- 
eral case taking into account the coupling between the 
lead explicitly, which include the information of all the 
cumulant of heat transfer. Also, the transient behavior 
and the long-time limit of CGF of heat transfer would 
be dealt with on an equal footing based on the NEGF 
method. 

This paper is organized as follows. We start in Sec. |TT] 
by introducing the model and the consistent quantum 
framework. And then we employ this framework to derive 
the CGF of heat transfer across the junction in Scc. lIIII In 
addition, we obtain a compact form of the long-time limit 
of the CGF in Sec. IIV1 based on which the steady-state 
fluctuation theorem (SSFT) is checked and the general- 
ized Caroli formula closely related to the first cumulant 
is recovered in Sec. |Vj At the end, we will briefly discuss 
Di Ventra's trick. 
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II. MODEL AND CONSISTENT QUANTUM 
FRAMEWORK 

We consider the lead-junction-lead model initially pre- 
pared in product state p ini = ^-X^ ® Tr( e -^c) 

— 7— Tr—jr^r- We can imagine that left lead (L). center 

junction (C), and right lead (R) in this model were in 
contact with three different heat baths at the inverse 
temperatures /3l = (^b^l) -1 , Pc = (fcsTc) -1 and 
Pr = {UbTr)~ 1 respectively, for time t < 0. At time 
t = 0, all the heat bath are removed, and coupling of 
the center junction with the leads Hlc = u^V LC uc 
and Hqr = UqV CR UR and the lead- lead coupling term 
Hlr = UjV Ur are switched on abruptly. Now the 
total Hamiltonian of the lead-junction-lead system will 
become 

Htot =Hl + He + Hr + Hlc + Hcr + Hlr, (1) 

where H a = ^p^,p a +^u^K a u a , a = L,C,R, represents 
coupled harmonic oscillators, and p a are col- 

umn vectors of transformed coordinates and correspond- 
ing conjugate momenta in region a. The superscript T 
stands for matrix transpose. 

In order to extract information on heat transfer, we 
introduce two-time observation protocol in the process 
of the evolution of the system; that is, at time t = 0~, 
we carry out the measurement of energy of the left lead 
associated with the operator Hl, obtaining the result 
to be the eigenvalue a of Hl, and measure it again at 
time t = tui, obtaining the eigenvalue b of Hl- Here the 
measurement is in the sense of quantum measurement 
of von Neumann^. This quantum history^, that is a 
sequence of quantum events at successive times, is repre- 
sented by the product projector Y = , where 
is a variant of the tensor product symbol 0, emphasizing 
that the factors in the quantum history refer to different 
times, and P^_is the projector on the energy eigenstatc 
of Hl with energy a measured at time t = 0~, similarly 
for Pt M . And the corresponding chain operator K (Y) is 
given by the expression K (Y) = Pf M U (tu,0~) Pq- in 
the case of a quantum history Y = Pq_ P^ M involving 
just two times, where U (4m>0 - ) is the time evolution 
operator of the full Hamiltonian H to t- It could be eas- 
ily verified that the joint probability distribution for the 
quantum history Y is 

Pr (Y) = (K*K) = Tr{p ini K*K} , (2) 

where the superscript f stands for transpose conjugate 
and from now on the angle brackets (...) simply denotes 
ensemble average with respect to p lnl . Here it is worth 
mentioning that this scheme for the joint probability dis- 
tribution could be readily extended to many-time mea- 
surement. 

Taking the commutator \p m% , Pq~] = into ac- 
count, we can verify that consistency conditions 



(K(Y),K(Y')) pini = Tr {p ini K (Y)* K(Y')} = for 

all F / F arc fulfilled, which guarantees that we are 
working with a consistent quantum framework. The 
consistent quantum framework, combining the axiom of 
probability with the Born rule of the quantum theory, is 
crucial to understand the probability aspect of the quan- 
tum theory^. Now we want to gain some insight into 
the specific consistent quantum framework we study from 
thermal transport point of view. 

As we will consider later, the generating function (GF) 
for heat transfer during the time Im is defined as 

2(0=E el?(Q ^ )pr ( P o a -©^). ( 3 ) 

a . b 

where the summation for a and b extends over all the 
eigenvalues of Hl- Then based on this definition, we 
simply take the derivative of Z (£) with respect to it; and 
then set £ = to obtain the average heat transfer out of 
left lead L 

Ql (*m) 

= E(a-6)Pr(P o a _0P ( b j/ ) (4) 
=Tr [p ini H L ] - Tr [p' n 'U (0", t M ) H L U (t M , 0")] , (5) 

where properties (Pq-) 2 = Pq- and $3 Pq_ = 1 f° r * nc 
projector , and similarly for P\ M are used. Also we 
use the key relation [p lm , Pq-] = 0, which assures that 
we are working with a consistent quantum framework. 
Thus, we immediately realize that 

dQ L {t M ) _ / dH L (t M ) \ (6) 
dt M \ dt M I ' 

of which the right-hand side is just the natural definition 
of thermal current II out of the left lead at time t = tu- 
Actually, if the initial density matrix does not commute 
with Pp_ , such as a steady-state density matrix, so that 
the framework we work with is not consistent, then the 
thermal current deriving from GF for the heat transfer 
will not be equal to the natural definition of it. 

III. CUMULANT GENERATING FUNCTION 

As was mentioned before, we proceed to study the GF 
for the heat transfer out of the left lead. According to 
the definition of the GF in Eq. we can obtain 

^(O =(Ut /2 (o-,t M ) u_ s/2 (t M ,o-)) , (7) 

where {tM, 0~) means evolution operator 

associated with the modified total Hamiltonian 
H~J /2 = e l ^/^ HL H tot e- i{ --^^ H ^, and similarly 
for (CP, ijv/)- Transforming to the interaction 

picture with respect to the free part of the modified 
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total Hamiltonian H = Hl + He + Hr, the GF for the 
heat transfer becomes 



2(0 = 

T e ~i So dru L (hx T +T)(V LC uc(r)+V LR u R (r))+uc(r)V CR u R (7 

(8) 



where T T is a r-ordering operator arranging operators 
with earliest r on the contour C (from - to <m and 
back to _ ) to the right, and a caret is put above oper- 
ators to denote their r dependence with respect to the 
free Hamiltonian such as uc (t) = e» ffcT !ice _ ^ Hc ' T , and 
x T = — £/2 with t = t + on the upper branch of the 
contour C, while x T = +£/2 with r = t~ on the lower 
branch. 

The key step to evaluate GF is to rewrite the exponent 
in Eg. (JSJ as -| J c dTidr 2 \u T (n) V (n, r 2 ) u (t 2 ), where 



(r) = [«l (&c T + t) w c (t) u_r (r)] , 








yLC 






V(n,T 2 ) = 







yCR 


5 (n , r 2 




yRL 








(9) 
(10) 



(r 1 ,r 2 ) = F(r 1 ,r 2 )) 



Here the generalized 8- function 8(ti,t 2 ) is simply coun- 
terpart of the ordinary Dirac delta function on the con- 
tour C, see, for example, Ref.— . 

Then expanding the exponential to perform a pertur- 
bation expansion and employing Fcynman diagrammatic 
technique, especially Wick's theorem and the linked clus- 
ter theorem, the CGF for the heat transfer can be ob- 
tained to be 



1 oo 

ln ^(o=^E Tr (^) 



n=l 
oo 



n=l 
1. 



-(Vg x y 

n 



~ {V re d9red) 



= - -Trln ( I - V red g 



red ) i 



(ii) 

(12) 
(13) 



and the notation Try r ) means trace both in real space 
index j and contour time r such as Tr( JjT ) [Vg x ] = 
J c dTi J c dT 2 TTj [V(Ti,T 2 )g x (t 2 ,ti)]; in the second 
equality, considering the structure of the last expression, 
instead of the full matrix we only need the hnitc reduced 
'ones, which make all the contribution to \nZ (£), and 



V red (ti,t 2 ) =8 (ri,r 2 ) 



V L 9 

red 



yLK 

red 



vh V™ 
vB v™ o 



(16) 



is the generalized 8- function 8{t\,t 2 ) times the reduced 
total coupling matrix obtained by deleting all the zero 
column and row vectors of the full one except for the 
possible zero vectors whose row or column indexes are 
the center (junction) ones, and g x ed is the corresponding 
submatrix of g x just like V red of V\ in the third equality 
a tilde above matrix means discretized contour-time ver- 
sion of the corresponding quantity such as [g x ed ] Tl T . n = 

\-9red\in ( T i^ T j)^ T j with an evenly spaced grid r, and Tj 
along the contour C and I is the identity matrix. 
Introducing the Dyson equation 

G(r,r') 

=9(t,t')+ / dn / dT 2 g(T,n)V (ti,t 2 ) G(t 2 ,t') 
Jc Jc 

(17) 



--9(t,t')+ I d n 
lc 



dr 2 G (r, n) V (n,r 2 ) g (r 2 , r') . 



(18) 



which actually defines G 



based on 



QLL qLC qLR' 
qCL qCC qCR 
qRL qRC qRR, 

g (t, t') = g x=a (t, t'), we easily realized that 

Gred — gred gredVredGred (1^) 
=3red + G re dV re dgred (20) 

still holds. And from now on, for notational simplicity, 
we omit all the subscript red with the understanding that 
these matrices are of finite dimensions in the real space 
domain. Thus, employing Egs. (|T5|) and ([2TJ)) along with 
the equality Trfnj4 = lndct A, we obtain 



where in the first equality 

g x (n,T 2 ) 

g L (hx Tl + T\ , hx T2 + t 2 ) 

g C (ri,T 2 ) 










g R (ri,r 2 ) 



(14) 



with the equilibrium contour-ordered Green functions 



9% (n,r 2 ) 



h r \Tr(e-^ 



a = L,C,R 



—T T [u?( ri )u% (r 2 )] 



(15) 



lnZ(0 = --lndet(J-d) 



d=[V LC V LR ] 



(21) 



~QCC QCR~ 




yCL- 


qRC qRR 




yRL 



~9t 



where g^ is the discretized contour-time version of 

st ( T i> T 2) = g L [hx Tl + n,hx T2 + r 2 ) - g L (n,T 2 ). 

If we introduce g~ x satisfying gJ^gL = I and employ 
Eqs. (QU) and ((20j), we can simplify the CGF for heat 
transfer further to be 
1 



In Z (£) = - | In dct [I - ~g- L l [G LL - g L j 

(22) 

which is valid in both transient and steady-state regimes. 
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IV. THE STEADY-STATE CGF 

Now, we proceed to evaluate the long-time limit of 
the CGF in Eg. (gU called the steady-state CGF. Trans- 
forming Eg. (|22[) from contour-time to real-time and then 
using the Keldysh rotation, which is essentially an or- 
thogonal transformation that 



A{t!,t 2 ) 



--0 T a z 



A(tf,t 2 ) A(tf,t 2 

A V'l i ^2 ) A [t 1 , t 2 

h) A K 



O 



A R (t u t 2 ) A a (t u t 2 ) 



(23) 
(24) 



with O 

1 

-1 



1 



1 1 

-1 1 



and the Pauli z matrix ov 



appearing due to the transition from contour- 



time to real-time, see, Eg. Rcf.— , we obtain 



where the notation TrQ t ) 
space j and real time t 
ft 



n V 



G 



LL 



9L 



(25) 



means trace both in real 



such as Tr 



(i,t) 



AB) = 



So d *i So ^Tr, \ A (*!, t a ) B (t 2 , h) . 

Before proceeding, a significant consideration comes 
that all kinds of real-time versions of the contour-time 
Green function G(t,t') defined in Eqs. (JHl) or ([IS]) 
are not necessarily time translationally invariant so that 
G LL (t, t') may not simply depend on the time difference 
t — t'. However, in the long-time limit, i.e., tm — > oo, 
the time translationally invariant part obtained from the 
lowest order of the Wigner transformation will dominate 
the CGFi£. It is equivalent to saying that G LL (t, if) = 
G LL (t — t') is time translationally invariant in the long- 
time limit (higher order terms of the product of the 
Wigner transformation have been ignored) . 

Consequently, setting t^i — > oo, and Fourier transform- 
ing Eq. ((23)) . we get 



1 



— In dct ( I - 

•no 27T 



where, 



9L 



9l 9l 
.0 9l 

.A 1 

9l = o 



G 



LL 



LL U LL 

G% L 



a — b a + b 
—a — b —a + b 



(26) 

(27) 
(28) 



g> (er^ -l),b = g< (e«"« - l) 



are all in frequency space. For NEGF notations and some 
general relations among Green's functions, we refer to 



To further simplify the steady-state CGF in Eq. (|26l) . 
/ A B N 



we use the formula det 



C D 



det {AD - BC) in case 



of [C, D] = to reduce the dimension of the matrix inside 
determinant by half. Therefore, the steady-state CGF is 
given by 

/oo 
— lndet{/-7^ M 
-oo 

x [(e^ - 1) f L (1 + f R ) + (e"^ - 1) f R (1 + f L )] }, 

(29) 



where Tg M =G r LR T R G a RL T L with 

-l 



{L,R} 



9{L,R} 



9{L,R} 



(30) 



is the transmission matrix and /{ l,r} = 

{exp (f3{L,R}hu) — l} 1 is the Bose-Einstein distri- 
bution function for phonons. 

In deriving it, we have used fluctuation dissipation the- 
orem e- 0LHw g2 M = gf [w] = f L (g r L - g a L ) along with 

Gll 

=G r LL (-</ifi) G a LL + G r LR (-if R t R ) G RL , (31) 

G r LL ~ G a LL ~ G r LL (gr 1 - gr 1 ) G a LL 
=G r LR (9 R - 1 -9 r R 1 )G RL , (32) 

due to (<?c) -1 — (gc) 1 = and the Langreth theorem^ 
acting on Eqs. (|19[) and (|2"0"1) . A computationally practi- 
cal closed equation for G r LR could be found in Ref^. The 
formulas Eq. (|29|) and Eq. (|22|) are our central results. 

Now we recover the classical version of the CGF for 
the heat transfer in harmonic networks without the lead- 
lead coupling, which was first derived in Ref£ using the 
Langevin equation method. To this end we simply set the 
lead-lead coupling V LR = in Eqs. (|T5J) and fl^b}, then 
use the Langreth theorem and the Fourier transforma- 
tion to obtain G r LR = g r L V LC G r CR = g r L V LC G r cc V CR g R 
along with G a RL = (G£ R ) f . 

After setting h — > and employing G RL [— w] = 
(G r LR H) T and f {L>R} [-u] = -f {LjR} [w] T , we can get 



hm 



J-TM 



/ — In dct <^ 

JO 27T I 

xk 2 B T L T R (iO [if+CSfl-jSi)]}, 



with 



T[w] =G 7 cc ri ? Gc C r L , 



(33) 

(34) 
(35) 



S {r,a} =y Ca ff {r,a} y aC 5 a = L R 
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V. THE SSFT AND CUMULANTS 

According to the steady-state CGF in Eq. (f29|). one 
could easily verify that the GC fluctuation symmetry^ 
Z (£) = Z (— £ + i (j3 R — (5 L )) is still satisfied in this gen- 
eral set-up with lead-lead coupling. And recall the def- 
inition of GF in Eq. ([3]), we know that the probability 
distribution for the heat transferred Ql is Pt(Ql) = 
^ JZc diZ (0 . Therefore, following the GC 

symmetry is the SSFT Pr (Ql) = e^ R ~^ QL Pr(-Q L ), 
which quantifies the violation of the second-law in the 
sense of probability. 

Also, the CGF can be used to evaluate cumulants. 
Here we only focus on steady-state cumulants of heat 
transfer. As illustrated in Sec. QJ the steady current is 
closed related to the first cumulant so that 

= J ^^-W Tr7 «H, (36) 

where, Z (0) = 1 is used. This generalized Caroli formula 
with lead-lead coupling was given very recently in Ref 
based on the definition of current directly, which gave 
us some valuable hints on the form of the steady-state 
CGF. The second cumulant describing the fluctuation 
of the heat transferred is obtained by taking the second 
derivative of steady-state CGF with respect to i£ and 
then setting £ = 0, which is 

((Ql)) =tM ^ ^ (M 2 { (h + f R + 2f L f R ) TrTg 

+ (/L-/*) 2 TrT g 2 }. (37) 

Higher-order cumulants are also systematically given by 
corresponding higher-order derivatives. 

After some experiences on first few order cumulants 
of heat transfer, we want to discuss the trick suggested 
by Di Vcntra mentioned in Sec. U that repartitioning the 
total Hamiltonian to avoid the inevitable coupling be- 
tween leads in real nanoscale or mesoscopic system when 
calculating steady current. Now whether this trick is 
applicable to the evaluation of higher-order cumulant 
(fluctuation) of the heat transfer in steady state boils 
down to checking whether TiTq old = TiTq new hold 
for all n less than or equal to the corresponding order 



of the cumulant one wants, where Tg.oU (Tg. new) is the 
transmission matrix before (after) repartitioning the to- 
tal Hamiltonian. Though giving a general verification is 
difficult, TyTq old = TtTq new ,Vn is indeed true in a one- 
dimensional central ring model, in which there is only one 
particle in the center junction connected with two semi- 
infinite spring chain leads and the interaction between 
the two nearest particles in the two leads respectively ex- 
ists (in this case, both Tc,oid and Tcnew are just a num- 
ber) . One step forward, if one think of CGF of heat trans- 
fer as the complete knowledge of the steady state, we can 
claim that the steady state is partition-independent after 
verification of TtTq^ = TiTq new , Vn or equivalently 
In Z id (£) = lnZ new (£) in Eq. (|29l) . Then we can partly 
answer a question raised by Caroli et al. regarding the 
(non)equivalence between the partitioned and partition- 
free approaches^, which recently was partly settled by 
explicitly constructing a non-equilibrium steady state 
through adiabaticaUy turning on an electrical bias be- 
tween the leads^i. 



VI. SUMMARY 

We examine the statistics of heat transfer during time 
tjvf in a general lead-junction-lead quantum system, in 
which coupling between leads has been taken into ac- 
count. To this end, a consistent quantum framework was 
introduced to derive the CGF valid in both transient and 
long-time regimes using the NEGF method. Also, the im- 
plication of consistency of the quantum framework was 
briefly discussed from thermal transport point of view. 
After that, a compact form of the steady-state CGF 
was obtained, following which the GC symmetry and the 
SSFT was verified. In addition, first few cumulants were 
given and generalized Caroli formula was recovered. Fur- 
thermore, some valuable hints with respect to the rig- 
orous proof for whether fluctuation of heat transfer in 
steady state is partition-independent have been offered. 
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